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We consider quantum superposition states in Bose-Einstein condensates. A de- 
coherence rate for the Schrodinger cat is calculated and shown to be a significant 
threat to this macroscopic quantum superposition of BEC's. An experimental sce- 
nario is outlined where the decoherence rate due to the thermal cloud is dramatically 
reduced thanks to trap engineering and "symmetrization" of the environment. We 
show that under the proposed scenario the Schrodinger cat belongs to an approximate 
decoherence-free pointer subspace. 
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CN ■ I. INTRODUCTION 

a '■ 

Microscopic quantum superpositions are an everyday physicist's experience. Macroscopic quantum superpositions, 
despite nearly a century of experimentation with quantum mechanics, are still encountered only very rarely. Fast 
decoherence of macroscopically distinct states is to be blamed In spite of that, recent years were a witness to 
an interesting quantum optics experiment g] on decoherence of a few photon superpositions. Moreover, matter- 
I ' wave interference in fullerene Cgo has been observed ||. Another experiment has succeeded in "engineering" the 
environment in the context of trapped ions Recently the first detection of a macroscopic Schrodinger cat state 
Q\ • in an rf-SQUID was reported J5j. All these successes tempt one to push similar investigations of basic quantum 
mechanics into the rapidly progressing field of Bose-Einstein condensation (BEC) of alkali metal atomic vapors |J 
The condensates can contain up to 10 7 atoms in the same quantum state. What is more, it is possible to prepare 
condensates in two different internal states of the atoms. Some of these pairs of internal levels are immiscible, and 
their condensates tend to phase separate into distinct domains with definite internal states Jt|. The immiscibility 
seems to be a prerequisite to prepare a quantum superposition in which all atoms are in one or the other internal 
state, — (\N, 0) + |0, N))/ \/2, where N is the total number of condensed atoms. There are at least two theoretical 
proposals how to prepare a macroscopic quantum superposition in this framework [||^]. Other proposals involve 
non-unitary evolution of the BEC towards the Schrodinger cat state by means of continuous quantum measurements 
I ■ |l(^[nj. Neither of these proposals addresses the crucial question of decoherence. 

This paper is a simplified and more pedagogical version of our previous work Jl2| . All technical details of the 
calculations, in particular the computation of the decoherence rate for the BEC superposition state, can be found in 
^ ' that reference. 



II. THE CONDENSATE 

Let us first make a short summary of some of the ideas for creating quantum superposition states of Bose-Eintein 
condensates. The methods that have been thus far proposed start with two weakly interacting dilute Bose condensates 
of atoms in different internal states. Let us call these internal states A and B. Atoms interact through s-wave 
collisions, characterized by a single parameter a which is known as the scattering length. We shall consider the case 
of repulsive interaction a > 0; it is known that for the opposite case the condensate is unstable above a critical 
total number of particles. It is assumed that the inter-scattering lengths for collisions A — A and B — B are the 
same, i.e. claa = &bb, which is in general different from the intra-scattering length ciab- When the self-energy of 
atom-atom interactions of the BEC is much less than the mode energy spacing, then one can treat the condensate in 
the two-mode approximation. Moreover, the condensate is shined with a appropriately chosen laser that introduces a 
Josephson-like coupling of strength A that interchanges internal atomic states in a coherent manner. The condensate 
two-mode Hamiltonian is 
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H c = e g (a)a + tfb) - \{a)b + tfa) + — (aVaa + tftfbb) + v c (a?tfab) , (1) 

Here e g is the energy of the lowest single-particle state in the trap (assumed to be the same for A and B), u c = 

4ir?i 2 aAA/'m and v c = inh dAB/m- 

This hamiltonian was studied in great detail in [^,^ 13 1. The precise quantum state of the condensate for all values 



of u c and v c was studied by means of a Wigner-like distribution on the two-mode Bloch sphere. It was found that 
when the repulsion between atoms in different states is bigger that the repulsion between atoms in the same state, 
that is when v c > u c , then the condensate tends to phase separate in space into distinct domains with definite internal 
states. This condition is known as the immiscibility condition. More importantly, the ground state of the system is a 
quantum superposition state. When the "purity" factor defined as e = (X/(v c — u c )N) N is much less than one, then 
the lowest energy subspace contains two macroscopic superpositions 

|+) = j=i[(J) N + (b f ) N ] 10, 0) = -±=(\N, 0) + |0, AO) , 

1 [( a t)*_( 6 t)"]| 0) 0) = i= (|AT,0)-|0,iV» . (2) 



The lower |+) and the higher |— ) states are separated by a small energy gap of E- — E + = N(u c — w c )elne. In the 
extreme case of A = (no laser applied) the two states are completely degenerate. If we just have e < 1, the |±) states 
in fact contain an admixture of intermediate states \N — 1, 1), . . . , |1, N — 1) such that their overlap is (+|— ) = e. For 
e ^> 1 they shrink to 

(at±&t)JV | 0>0) . (3) 

which does not correspond to a macroscopic superposition. From now on we assume the pure case, e<l. 

The proposal of Gordon and Savage H for obtaining these superposition states is similar to a typical experiment in 
quantum optics with non-linear systems. In the BEC case, the non-linearity is provided by the collisional interactions 
and the Josephson coupling. The idea is to prepare an initial state with all the atoms in the same internal state and 
then turn on the Josephson coupling for some amount of time adequately chosen. After the Josephson coupling is 
turned off, the quantum state of the system has been modified, and a Schrodinger cat state has been formed. Other 
proposals for obtaining basically the same kind of superposition use a continuous quantum measurement process [fll| . 

III. THE THERMAL CLOUD 

As we have already mentioned, neither of the proposals for generating quantum superpositions in BEC addresses 
the crucial question of decoherence. In fact, the condensate is an open quantum system which is in contact with an 
environment of non-condensed thermal particles. This interaction may be responsible for the loss of coherence between 
the components of the quantum superposition states of the previous section. If the decoherence time is very small, 
then the existence of these states in BEC would be merely of academic interest, since the expectations of observing 
them in the lab would vanish. Therefore it is important to understand how the thermal cloud affects the longevity of 
the BEC cats. 

The Hamiltonian for the dilute environment formed by non-condensed particles is 

H E =J2bs(ala s + blb s )-X(alb s + bta s )] . (4) 

s 

e s is the single particle energy of level s. After the transformation 

s + b s _ As 

V2 ' S ~ V2 



G a s + b s a s ~b s 

£>s = t=— , O s = ■=— (5) 



adopts a diagonal form 

H E = Y / [(ts-X)StS s + (e s + \)OtO s ] . (6) 

s 

In Figure 1 we make a schematic plot of the energy levels of the non-condensed system before (A = 0) and after 
(A 7^ 0) applying the Josephson coupling. For simplicity we consider an isotropic harmonic trap with a dip in its 
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center, where the condensate particles will be localized (we shall return to this point later). We see that when the 
laser coupling is applied, the states S s (which are symmetric under the exchange a <-> b) and the states O s (which 
are antisymmetric under the same exchange) feel the same, but shifted, trapping potentials. The two ladders of 
energy eigenstates of the environment are separated by a gap equal to 2A. Each of the levels is occupied according 
to the thermal Bose distribution. When the gap is much bigger that the typical thermal energy 2A 3> k^T, then 
the antisymmetric states are unoccupied. We will see in the next section that this is a prerequisite for reducing 
decoherence. 



Condensate 




Non 
Condensate 



FIG. 1. Environmental single-particle energy levels before and after turning on the Josephson-like coupling. 



IV. CONDENSATE-NON CONDENSATE INTERACTIONS 



Condensed atoms interact with non-condensed atoms via two-body collisions. The exact interaction hamiltonian V 
can be easily obtained from the total (Gross-Pitaevskii like) hamiltonian 



H 



(7) 



v^aMUb) + + V^V^ + U{r)^ A <f> A - WUb] + {A <-> B} 

after splitting the quantum fields 4>a and <\>b into condensate and bath parts, i.e. 

4>a{x) = a g(x) + ^ a s u s (x) , 

S 

4> B {x) = b g(x) + ^2 b sU s (x), (8) 



where g(x) is the wave function of the condensate localized in the dip. After introducing the symmetric and anti- 
symmetric operators S s and O s for the thermal bath, the interaction hamiltonian can be splitted into two terms, 
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V = Vs + Vo, which are respectively symmetric and antisymmetric under the interchange a «-> b of system operators. 
For the sake of conciseness, we shall refrain from writing down the whole expressions for these two interaction hamil- 
tonians. Suffice it to say that each of them has two possible energy-conserving interaction verteces, both proportional 
to the interaction couplings norw: 

• There are inelastic two-body collisions that do not conserve the number of condensed particles (condensate 
feeding and depletion processes). There are in turn two possible diagrams for these processes, and they are 
shown in Figure 2. The first one is proportional to the occupation number of initial non-condensed states; since 
each occupation number is proportional to fugacity z = (where (5 is the inverse temperature of the thermal 
bath and /x its chemical potential), then the diagram is 0(z 2 ). The second one is proportional to the occupation 
number of the initial non-condensed state (that is, proportional to z), and energy conservation brings about 
another fugacity factor, so the final result is that this diagram is also 0(z 2 ) |I4| , ^5| . 




FIG. 2. Inelastic two-body collisions between condensate and non-condensate particles 

• There are elastic two-body collisions that conserve the number of condensed particles, and are shown in Figure 
3. These diagrams are proportional to fugacity, O(z) p^ , [l5| ]. 




FIG. 3. Elsatic two-body collisions between condensate and non-condensate particles 

When fugacity is small (which occurs when the gap between the condensate and non-condensate single-particle 
levels is bigger that the thermal energy k^T), then the elastic processes dominate over the inelastic ones. The elastic 
collisions are also the most relevant ones in the computation of the decoherence rate, since we shall see that their 
contribution scales as TV 2 (with N the total number of particles in the condensate) , whereas the inelastic contribution 
scales as N. For these reasons, from now on we shall concentrate only on the elastic processes. 

V. DECOHERENCE-FREE POINTER SUBSPACE IN BEC 

As we have mentioned in Section III, when 2A 3> k^T the antisymmetric environmental states are nearly empty. 
Only the symmetric states are occupied. Since these latter states cannot distinguish between A and B, all collisions 
that involve symmetric non-condensed states will not destroy the quantum phase coherence between the Schrodinger 
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cat's components. We refer to this condition as the perfect symmetrization limit. In this case, the states |±) = 
( | iV, 0) ± 1 , A^) ) / \/2 span a decoherence-free pointer subspace (also known in the literature as decoherence-free subspace 
or, for short, DFS) of the Hilbert space, since they have degenerate eigenvalues of the interaction hamiltonian V [p~6| — p~8|| . 
Any state of that subspace can be written as a\N, 0) + /3|0, N), with a and (3 complex numbers. If / P[ a |7V,o)+/3|o,./v)] 
denotes a projector onto that subspace, then 

[V,V[ a \N tO )+0\o,N)]] = 0, (9) 

which means that any quantum superposition a\N, 0) + /3|0, N) in the subspace is an eigenstate of the interaction 
hamiltonian (a perfect pointer state), and as such will retain its phase coherence and last forever. The interaction 
Hamiltonian between the condensate and the thermal cloud is a sum of products of condensate operators and environ- 
mental operators. Only terms with symmetric environmental operators are relevant because the antisymmetric states 
are empty. The total Hamiltonian is symmetric with respect to A «-* B so, to preserve this symmetry, the relevant 
terms with symmetric environmental operators also contain symmetric condensate operators. The argument simplifies 
a lot for small fugacity where there is only one leading term with the Na + Nb condensate operator. The states |±) 
are its eigenstates with the same eigenvalue N. They are also (almost) degenerate eigenstates of the condensate 
Hamiltonian build out of Na.b- The coherent transitions A <-> B break this degeneracy of |±) but the difference 
of their eigenfrequencies is negligible as compared to the usual condensate lifetime of ~ 10s. In the next-to-leading 
order in fugacity there are symmetric interaction terms which change the number of condensed atoms. These terms 
drive the \N, 0) and |0, AT) states into slightly "squeezed-like" states \S, 0) and |0,S) respectively p^ |. There are also 
terms which exchange A with B. They give each state a small admixture of the opposite component. Superpositions 
of these are still decoherence-free pointer subspaces - there are no relevant antisymmetric operators to destroy their 
quantum coherence. 



VI. DECOHERENCE COMES INTO THE SCENE 

When the antisymmetric environmental states begin to be occupied, then the commutation relation (||) is only 
approximate and states within the subspace will decohere. Hence we must face the problem of calculating decoherence 
rates due to interactions with the environment. 

We shall first consider the case of A = 0, which corresponds to no Josephson-coupling being applied, and for which 
the two states |±) are exactly degenerate with respect to the condensate hamiltonian. In this case there is no difference 
between symmetric and antisymmetric environmental states. After some long but straightforward calculation which 
involves computing the master equation for the reduced density matrix of the condensed particles (see ]l^] for details), 
we find a lower bound for the decoherence rate 

3 (a 2 N E 



*dec > 16t j v T N A , (10) 



where A^ is a number of condensed atoms, vt — \j2k^,T '/m is a thermal velocity in the noncondensed thermal cloud, 
a is a scattering length, V is a volume of the trap, and Ae is a number of atoms in the thermal cloud, 

N* * e& ( . (11) 



This is a lower estimate since we have only considered terms to leading order in fugacity O(z) and to leading order 



in condensate size 0(N 2 ). Next-to-leading order terms are 0(z 2 ) and 0(N), in agreement with Refs. |14 15[ ], so they 
were neglected here. Even without going into details of our derivation it is easy to understand where a formula like 
Eq.(|l0|) comes from. A^ 2 is the main factor which makes the decoherence rate large. It comes from the master equation 
of the Bloch-Lindblad form p ~ — [N\ — N&, \N\ — A^b,p]], with A and B the two internal states of the atoms. A^ 2 
is the distance squared between macroscopically different components of the superposition (|AT, 0) + [0, N))/ \/2 - the 
common wisdom reason why macroscopic objects are classical [Q. The factor in brackets in Eq.([l0|) is a scattering 
rate of a condensate atom on noncondensate atoms - the very process by which the thermal cloud environment learns 
the quantum state of the system. 

Let us estimate the decoherence rate for a set of typical parameters: T — 1/xK, lj — 50Hz, and a = 3 . . . 5nm. 
The thermal velocity is «t ~ 10 _2 m/s. The volume of the trap can be approximated by V — 47ra^ cturn /3, where 
flrcturn = \J 2k-gT / mid 2 is b, return point in a harmonic trap at the energy of k^T . We estimate the decoherence time 
as t dec « 10 5 sec/(AT E Af 2 ). For N E = 10° . . . 10 4 and N = 10... 10 7 it can range from 1000s down to 10" 13 s. For 
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TV = 10 our (over-)estimate for idee is large. However, already for N = 1000 and = 10 (which are still within 
the limits of validity of the two- mode approximation) tdec shrinks down to milliseconds. Given that our i^cc is an 
upper estimate and that big condensates are more interesting as Schrodinger cats, it is clear that for the sake of cat's 
longevity, one must go beyond the standard harmonic trap setting. 



VII. TRAP ENGINEERING 



From Eqs.( |io| , pd| ) it is obvious that the decoherence rate depends a lot on temperature and on chemical potential. 
The two factors strongly influence both N-e and vt- Both can be improved by the following scenario, which is 
a combination of present day experimental techniques. In the experiment of Ref. |l9| a narrow optical dip was 
superposed at the bottom of a wide magnetic trap. The parameters of the dip were tuned so that it had just one 
bound state. The gap between this single condensate mode and the first excited state was 1.5/iK, which at T = 1/iK 
gives a fugacity of z = exp(— 1.5). We need the gap so that we can use the single mode approximation. At low 
temperatures, the gap results in a small fugacity, which is convenient for calculations. We propose to prepare a 
condensate inside a similar combination of a wide magnetic and a narrow optical trap (or more generally: a wide 
well plus a narrow dip with a single bound mode) and then to open the magnetic trap and let the noncondensed 
atoms disperse. The aim is to get rid of the thermal cloud as much as possible. A similar technique was used in the 
experiment of Ref. ^(J (see Figure 4) . 




FIG. 4. First step in engineering the trap: an optical dip is superimposed on a big magnetic trap, which is later opened so 
that most of the thermal atoms disperse away. 



Let us estimate the ultimate limit for the efficiency of this technique. At the typical initial temperature of 1/iK the 
thermal velocity of atoms is 10 _2 m/s. An atom with this velocity can cross a 1/zm dip in 10~ 4 s. If we wait for, say, Is 
after opening the wide trap, then all atoms with velocities above 10 _6 m/s will disperse away from the dip. A thermal 
velocity of 10~ 6 m/s corresponds to the temperature of 10 _8 ^K. As the factor Nevt ~ T 7 / 2 in Eq.(|i~0|), then already 
Is after opening the wide trap the decoherence rate due to non-condensed atoms is reduced by a factor of 10~ 28 ! 

It is not realistic to expect such a "cosmological" reduction factor. The "dip" which is left after the wide harmonic 
trap is gone could be, for example, a superposition of an ideal dip plus a wide shallow well (which was a negligible 
perturbation in presence of the wide harmonic trap). The well would have a band of width AE of bound states which 
would not disperse but preserve their occupation numbers from before the opening of the wide trap. They would stay 
in contact with the condensate and continue to "monitor" its quantum state. Even if such a truncated environment 
happens to be already relatively harmless, there are means to do better than that. 

Further reduction of the decoherence rate can be achieved by "symmetrization" of the environmental states. To 
this end we propose to turn on the Josephson-like coupling A in order to induce coherent transitions from states A and 
B. This has already been achieved experimentaly [pl| ; in their case A « 1kHz. When AE <C 2A wc arc in the perfect 
symmetrization limit, and hence no decoherence takes place. Indeed, the symmetric and antisymmetric AS-bands 
of states can be visualized as two ladders shifted with respect to each other by 2A. In other words, the two sets of 
states feel the same, but shifted, trapping potentials. When AE <C 2A, then the antisymmetric O s 's are nearly empty 
since they can evaporate into symmetric states and then leave the trap . The symmetric S s 's cannot distinguish 
between A and B so they do not destroy the quantum coherence between the Schrodinger cat's components. After 
symmetrization, N-e in Eq.(|i"o|) has to be replaced by the final number of atoms in the antisymmetric states only; 

atO „ J n x exp((/i - X)/k B T), for 2A < AE 

iVE ~ \ 0, for 2A > AE (12 > 

Here T is the temperature before opening the wide trap. n\ is the number of antisymmetric bound states which 
remain within the AS-band of symmetric states. Atoms in these antisymmetric bound states cannot disperse away. 
For 2A > AE this number n\ is zero and there is no decoherence from the thermal cloud. 
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In Figure 5 we make a schematic graph of the steps we propose to engineer the trap in order to reduce decoherence 
due to the thermal cloud. 




FIG. 5. Second step in engineering the trap: After the magnetic trap is opened, there is still a band AE of thermal states in 
the "mouth" of the dip. Turning on the Josephson coupling A shifts the energy levels of symmetric and antisymmetric thermal 
states, and the latter become empty. 

VIII. OTHER SOURCES OF DECOHERENCE 

Interaction with the thermal atoms is not the only source of decoherence for the quantum superposition state of 
the condensate. Among other possible sources we can mention: 

• Ambient magnetic fields: the condensed atoms have magnetic moments. If the magnetic moments of A and B 
were different the magnetic field would distinguish between them and would introduce an unknown phase into 
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the quantum superposition, thus rendering its underlying coherent nature undetectable. Fortunately the much 
used \F, mp) — |2, 1), |1, —1) states of 87 Rb have the same magnetic moments. For them the magnetic field is a 
"symmetric" environment. 

• Different scattering lengths: there is a typical « 1% difference between the A — A and B — B scattering lengths. 
The condensate Hamiltonian is not perfectly symmetric under A <-> B. Even for a perfectly symmetrized 
environment, symmetric environmental operators couple to not fully symmetric condensate operators. This 
means that for the 1% difference of scattering lengths symmetrization can improve decoherence time by at most 
two orders of magnitude as compared to the unsymmetrized environment. 

• Three-body losses: due to collisions involving three particles, the condensate is not stable but loses atoms via 
recombination into molecular states. Hence a condensate has a finite lifetime, of the order of 10 . . . 20s. The 
atoms which escape from the condensate carry information about its quantum state. They destroy its quantum 
coherence. In the experiment of Ref. [|19| the measured loss rate per atom was 4/s for ./V = 10 7 or around 1 atom 
per 10~ 7 s. The last rate scales like so already for N — 10 4 just one atom is lost per second; decoherence 
time is Is. One possibility for increasing this decoherence time is to increase slightly the dip radius. The loss 
rate scales like density squared so an increase in the dip width by a factor of 2 reduces the loss rate by a factor 
of 2 6 = 64. 



IX. DISCUSSION 

The aim of this paper was to discuss the "longevity" of Schrodinger cats in BEC's. We have shown that while in the 
standard traps decoherence rates are significant enough to prevent long-lived macroscopic superpositions of internal 
states of the condensate, the strategy of trap engineering and symmetrization of the environment will be able to deal 
with that issue. 

What remains to be considered is how one can generate such macroscopic quantum superposition, and how one can 
detect it. The issue of generation was already touched upon in Refs. We have little to add to this. However, 

in the context of the Gordon and Savage proposal, it is fairly clear that the time needed to generate the cat state 
would have to be short compared to the decoherence time. If our estimates of Eq.(|l^) are correct, symmetrization 
procedure appears necessary for the success of such schemes. 

Detection of Schrodinger cat states is perhaps a more challenging subject. In principle, states of the form (\N, 0) + 
\0,N))/>/2 have a character of GHZ states, and one could envision performing measurements analogous to those 
suggested in J27J and carried out in j28|, where a 4-atom entangled state was studied. However, this sort of parity- 
check strategy, appropriate for TV < 10, is likely to fail when N is larger, or when (as would be the case for the 
"quasi-squeezed" states anticipated here [p2[ ) N is not even well defined. 

A strong circumstantial evidence can be nevertheless obtained from two measurements. The first one would consist 
of a preparation of the cat state, and of a measurement of the internal states of the atoms. It is expected that in 
each instance all (to within the experimental error) would turn up to be in either A or B states. However, averaged 
over many runs, the number of either of these two alternatives would be approximately equal. Decoherence in which 
the environment also "monitors" the internal state of the atoms in the A versus B basis would not influence this 
prediction. We need to check separately whether the cat state was indeed coherent. To do this, one could evolve 
the system "backwards". However, this is not really necessary. For, as Gordon and Savage point out, when, in their 
scheme, we let the system evolve unitarily for more or less twice the time needed for the generation of the cat state, 
it will approximately return to the initial configuration. Thus, we can acquire strong evidence of the coherence of the 
cat provided that this unitary return to the initial configuration can be experimentally confirmed. 

These are admittedly rather vague ideas, which serve more as a "proof of principle" rather than as a blueprint for 
an experiment. Nevertheless, they may, we hope, encourage more concrete investigation of such issues with a specific 
experiment in mind. 
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